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Abstract 

The  McCallum  Projection  for  Cylindrical  Algebraic  Decomposition 
(CAD)  produces  a  smaller  projection  factor  set  than  previous  projections, 
however  it  does  not  always  produce  a  sign-invariant  CAD  for  the  set  of 
input  polynomials.  Problems  may  arise  when  a  (fc  -I-  l)-level  projection 
factor  vanishes  identically  over  a  fc-level  cell.  According  to  McCallum’s 
paper,  when  this  happens  (and  k+1  is  not  the  highest  level  in  the  CAD)  we 
do  not  know  whether  the  projection  is  valid,  i.e.  whether  or  not  a  sign- 
invariant  CAD  for  the  set  of  input  polynomials  will  be  produced  when 
lifting  is  performed  in  the  usual  way.  When  the  fc-level  cell  in  question 
has  dimension  0,  McCallum  suggests  a  modification  of  the  lifting  method 
that  will  ensure  the  validity  of  his  projection,  although  to  my  knowledge 
this  has  never  been  implemented. 

In  this  paper  we  give  easily  computable  criteria  that  often  allow  us 
to  conclude  that  McCallum’s  projection  is  valid  even  though  a  projection 
factor  vanishes  identically  over  a  cell.  We  also  improve  on  McCallum’s 
modified  lifting  method. 

We’ve  incorporated  the  ideas  contained  in  this  paper  into  QEPCAD, 
the  most  complete  implementation  of  CAD.  When  McCallum’s  projec¬ 
tion  is  invalid  because  of  a  projection  factor  not  being  order-invariant 
over  a  region  on  which  it  vanishes  identically,  at  least  a  warning  message 

‘This  article  was  originally  posted  on  my  website  in  September  2001  as  MQTS2001 . 1  .ps  .gz. 

It  is  now  being  reformatted  and  republished  as  a  USNA  CS  Department  Technical  Report. 

All  references  to  “current”  are  to  be  taken  as  “current  as  of  1  September,  2001. 
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ought  to  be  issued.  Currently,  QEPCAD  may  print  warning  messages 
that  are  not  needed,  and  may  fail  to  print  warning  messages  when  they 
are  needed.  Our  implementation  in  QEPCAD  ensures  that  warning  mes¬ 
sages  are  printed  when  needed,  and  reduces  the  number  of  times  warning 
messages  are  printed  when  not  needed.  Neither  McCallum’s  modified 
lifting  method  nor  our  improvement  of  it  have  been  implemented  in  QEP¬ 
CAD  —  the  design  of  the  system  would  make  implementing  such  a  feature 
quite  difficult. 


1  Introduction 

The  McCallum  Projection  [11]  (and  the  Improved  McCallum  Projection  [3])  rep¬ 
resents  a  huge  improvement  over  the  original  projection  [6]  for  CAD  construc¬ 
tion,  as  the  projection  factor  set  it  produces  is  much  smaller.  Unfortunately, 
it  also  involves  a  slightly  more  complicated  idea  of  lifting.  In  particular,  lifting 
becomes  difficult  when  a  projection  factor  vanishes  identically  over  some  set  (for 
example,  as  {x  +  y)z  —  (x  —  y)  vanishes  identically  over  (x  =  0,y  =  0)).  Accord¬ 
ing  to  McCallum’s  original  paper,  when  a  projection  factor  vanishes  identically 
over  a  region  of  dimension  greater  than  zero,  the  projection  fails  (unless  that 
projection  factor  is  of  the  highest  level).  If  a  projection  factor  vanishes  identi¬ 
cally  over  a  region  of  dimension  zero,  then  a  “delineating  polynomial”  must  be 
added  to  the  lift  basis  in  that  step  (once  again,  unless  that  projection  factor  is 
of  the  highest  level).  The  key  idea  is  that  iterating  the  lifting  process  requires 
not  simply  sign-invariance,  but  the  stronger  property  of  order-invariance,  and 
while  a  polynomial  is  sign-invariant  over  a  region  over  which  it  is  identically 
zero,  it  is  not  necessarily  order-invariant. 

In  this  paper,  we  describe  more  precise  (and  yet  still  easily  computable)  criteria 
for  determining  when  order-invariance  is  required,  when  a  delineating  polyno¬ 
mial  needs  to  be  added  to  ensure  order-invariance,  and  what  that  delineating 
polynomial  should  be.  These  improvements  are  of  considerable  practical  im¬ 
portance.  When  the  McCallum  projection  fails,  either  direct  and  substantial 
work  by  hand  is  required  of  the  user,  or  the  McCallum  Projection  needs  to  be 
abandoned  in  favor  of  far  costlier  projections. 

We  consider  in  this  paper  the  QEPCAD  system,  which  provides  the  most  com¬ 
plete  implementation  of  CAD.  When  the  McCallum  Projection  requires  the 
addition  of  a  delineating  polynomial,  QEPCAD  fails!  ^  Not  only  doesn’t  QEP¬ 
CAD  currently  include  facilities  for  adding  delineating  polynomials,  but  its  ar¬ 
chitecture  makes  it  extremely  difficult  to  add  such  a  feature.  Thus,  apart  from 

^In  fact,  the  current  situation  is  worse:  QEPCAD  warns  the  user  if  a  projection  factor  not 
of  the  highest  level  vanishes  identically  over  a  positive  dimensional  cell,  but  simply  ignores 
cases  in  which  a  projection  factor  vanishes  over  a  0-dimensional  cell!  It  may  be  possible  that 
QEPCAD  gives  incorrect  answers  as  a  result. 
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any  performance  benefit  of  not  adding  unnecessary  delineating  polynomials, 
identifying  only  those  cases  in  which  they  are  truly  necessary  is  critical. 

These  improvements  have  been  implemented  in  QEPCAD,  with  the  result  that 
for  the  first  time  the  user  is  able  to  use  the  McCallum  Projection  and  be  assured 
that,  in  the  absence  of  error  messages  (such  as  “A  Delineating  Polynomial  Must 
Be  Added”),  the  McCallum  Projection  is  definitely  valid.  When  a  projection 
factor  P  vanishes  identically  over  cell  c: 

1.  We  determine  whether  order-invariance  is  really  required  for  P  (described 
in  Section  4),  if  not  QEPCAD  simply  continues. 

2.  When  c  is  0-dimensional  we  determine  whether  a  delineating  polynomial 
must  be  added  to  ensure  the  order-invariance  of  P  over  c  (described  in 
Section  2),  if  so,  we  print  an  error  message,  otherwise  QEPCAD  simply 
continues. 

3.  When  c  has  positive  dimension,  we  are  sometimes  able  to  determine  that 
no  delineating  polynomial  is  needed  to  ensure  the  order-invariance  of  P 
over  c  (described  in  Section  3)  and  in  this  case  QEPCAD  simply  continues 
—  otherwise  we  print  an  error  message. 


This  paper  assumes  that  the  reader  is  familiar  with  the  McCallum  projection 
and  the  basic  methods  of  CAD  and  quantifier  elimination  by  CAD  as  explained 
in,  for  example,  [7]  or  [2]. 


2  The  zero-dimensional  case 


Suppose  that  P{xi, . . . ,  Xk+i)  is  a  projection  factor  that  vanishes  identically  at  a 
point  a  =  (oi, . . . ,  Ofc),  i.e.  P{ai, . . . ,  Xk+i)  =  0.  In  this  section  we  consider 
the  problem  of  constructing  a  decomposition  of  the  line  a  x  K  into  regions  in 
which  the  order  of  P  is  constant. 


2.1  Minimal  delineating  polynomials 

When  projection  factor  P{xi , . . . ,  Xk+i)  vanishes  identically  over  a  point  a  G  M*, 
we  know  that  P  has  order  greater  than  0  in  a  x  M,  but  we  don’t  know  whether 
or  not  the  order  of  P  is  constant  in  that  line.  Suppose  t  is  the  smallest  index 
for  which  at  least  one  of  the  t-order  partials  of  P  is  not  identically  zero  over  a, 
and  let  Dp  be  the  set  of  all  t-order  partials  of  P.  We  know  that  the  order  of 
P  in  a  X  M  is  t  almost  everywhere,  but  at  finitely  many  points  the  order  may 
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be  greater  than  t.  In  fact,  the  order  of  P  in  a  x  M  is  greater  than  t  at  exactly 
those  points  at  which  all  the  elements  of  Dp  are  zero. 

We  would  like  to  introduce  a  “delineating  polynomial”  to  define  our  decompo¬ 
sition  of  a  X  M.  A  delineating  polynomial  is  a  nonzero  polynomial  in  IR[a;fc+i] 
whose  roots  include  the  a;fc+i-coordinate  of  every  point  in  a  x  M  at  which  the 
order  of  P  is  greater  than  t.  McCallum  points  out  that  any  element  of  Pp  that 
is  not  identically  zero  over  a  may  be  used  as  a  delineating  polynomial.  How¬ 
ever,  while  the  vanishing  of  such  a  polynomial  is  a  necessary  condition  for  an 
increase  in  the  order  of  P,  it  is  not  a  sujficient  condition.  All  elements  oi  Dp 
must  vanish  at  a  point  in  a  x  M  for  the  order  to  increase.  A  polynomial  which 
vanishes  at  exactly  the  the  x^+i-coordinate  of  every  point  in  a  x  M  at  which 
the  order  of  P  is  greater  than  t  will  be  called  a  minimal  delineating  polynomial. 
The  vanishing  of  a  minimal  delineating  polynomial  provides  a  necessary  and 
sufficient  condition  for  the  increase  in  order  of  P. 

Constructing  a  minimal  delineating  polynomial  is  quite  easy.  If  we  define 
S  =  {Q{ai, . . . ,  ttfc,  Xk+i)  I  Q  G  Dp} 

the  order  of  P  is  greater  than  t  at  exactly  the  points  of  the  form  (oi, . . . ,  /?), 

where  /3  is  a  zero  of  all  the  elements  of  S'.  If  G  is  the  GCD  of  the  elements  of  S, 
then  the  roots  of  G  are  exactly  the  common  roots  of  the  elements  of  S  —  i.e. 
G  is  a  minimal  delineating  polynomial! 

Thus,  instead  of  choosing  one  of  the  non-zero  elements  of  Pp  and  using  it  as  a 
delineating  polynomial  (as  suggested  in  [11]),  we  use  the  GCD  of  all  non-zero 
elements  of  Dp  as  a  delineating  polynomial. 


2.2  In  practice 

One  of  the  practical  benefits  of  computing  a  minimal  delineating  polynomial 
when  faced  with  a  projection  factor  P(a:i, . . . ,  Xfc+i)  that  vanishes  identically 
over  a  point  a  G  is  that  one  often  finds  that  the  minimal  delineating  poly¬ 
nomial  is  constant,  meaning  that  no  polynomial  needs  to  be  added  at  all!  Some 
examples  illustrate  this  point. 

Example  1  Consider  the  polynomial  P  =  {x  +  y  —  l)z  +  {y  —  1)  and  the 
point  a  =  (0, 1).  P  vanishes  identically  over  a,  so  according  to  the  McCallum 
projection  we  need  to  add  a  delineating  polynomial  to  decompose  a  x  M  into 
regions  in  which  P  is  order-invariant. 

Dp  =  {Px  =  z,Py  =  z+l,P:,  =  x  +  y  -1} 

Thus,  S  =  {z,z+l,  0}  and  G,  the  GCD  of  the  elements  of  S,  is  1.  This  means 
that  even  though  P  vanishes  identically  over  a,  the  order  if  P  is  invariant  in 
a  X  M.  Thus,  no  delineating  polynomial  is  needed! 
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Example  2  In  the  “X  Axis  Ellipse”  problem  (see  Section  6.1),  a  well-known 
problem  in  the  literature,  the  polynomial 

P{a,  b,  c,  x)  =  -  a^x^  -  2b^cx  +  b^c^  -  a^b^  +  c? 

appears  as  a  projection  factor  (it  is  the  resultant  of  two  initial  polynomials). 
The  point  a  =  (1,1,0)  is  a  zero-dimensional  cell  in  the  CAD  constructed  for 
this  problem,  and  P  vanishes  identically  over  a. 

Pa  =  —2ax^  —  2ab^  +  2a,  Pf,  =  2bx^  —  4bcx  +  2bc^  —  2ba^, 

Pc  =  —2b^x  +  26^c,  Px  =  2b^x  —  2a?x  —  26^c 

Thus,  S  =  {—2x^,2x^  —  2,— 2a:,  0}  and  G,  the  GCD  of  the  elements  of  S,  is 
l.Thus,  no  delineating  polynomial  is  needed! 


2.3  Must  non-constant  minimal  delineating  polynomials 
always  be  added? 

From  the  perspective  of  QEPCAD,  if  a  delineating  polynomial  really  does  need 
to  be  added,  we  are  out  of  luck.  Moreover,  a  considerable  amount  of  work  will  be 
required  to  modify  the  program  to  allow  the  addition  of  delineating  polynomials. 
So,  is  all  hope  lost  when  we  compute  a  minimal  delineating  polynomial  and  it 
is  non-constant?  Not  always! 

Let  P{xi, . . .  ,Xk+i)  be  our  familiar  projection  factor  that  vanishes  identically 
over  a  point  a  G  When  we  lift  over  a  during  CAD  construction,  P  is 

typically  one  of  many  {k  +  l)-level  projection  factors,  and  we  compute  a  lift 
basis,  namely  a  squarefree  basis  constructed  from  substituting  a  into  all  (fc-l- 1)- 
level  projection  factors,  which  defines  our  decomposition  of  a  x  M.  If  this  lift 
basis  already  “contains”  the  minimal  delineating  polynomial  we’ve  computed 
(i.e.  the  squarefree  part  of  the  minimal  delineating  polynomial  divides  the 
product  of  the  basis  polynomials),  then  there’s  nothing  to  “add”! 

This  might  seem  unlikely,  but  in  fact  it’s  not.  In  Anai’s  Problem  (see  Sec¬ 
tion  6.2),  we  have  a  projection  factor  P{s,  z,  xi,X2)  =  SX2  —  xf.  The  point  a  = 
(0,-1,  0)  is  a  zero-dimensional  cell  in  the  CAD  constructed  for  this  problem,  and 
P  vanishes  identically  over  a.  Dp  =  {Pg  =  X2,Pz  =  0,Pa;i  =  —2xi,Px2  =  s}) 
and  thus  S  =  {x2,  0,0,0}  and  the  minimal  delineating  polynomial  is  X2-  From 
this  we  see  that  the  order  of  P  in  a  x  M  is  1  everywhere  but  (0,  —1,  0, 0),  where 
it  is  higher  (order  2,  specifically).  It  seems  that  we  must  add  a  delineating 
polynomial.  However,  X2  is  itself  a  projection  factor,  and  thus  a  part  of  the  lift 
basis.  Therefore,  the  minimal  delineating  polynomial  doesn’t  need  to  be  added, 
it’s  already  there!  [Note  that  as  a  quantifier  elimination  problem,  the  vanishing 
projection  factor  would  not  cause  difficulties  because  it  is  at  the  highest  level, 
so  sign-invariance  suffices.  What  the  above  argument  shows  is  that  what  we 
actually  get  an  order-invariant  decomposition.] 
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3  Dimension  greater  than  zero 


In  Section  2.1  we  deal  with  the  vanishing  of  projection  factor  P  over  a  point  a.  by 
considering  the  system  of  all  t-order  partials  of  P.  We  often  find  that  the  system 
is  inconsistent,  and  therefore  determine  that  P  is  actually  order-invariant  over 
a.  In  this  section  we  try  to  apply  the  same  ideas  to  the  situation  in  which  P 
vanishes  over  a  region  of  higher  dimension. 


3.1  In  theory 

If  c  is  a  fc- level  cell  in  a  CAD,  it  is  a  connected,  open  subset  of  some  variety 
defined  by  all  the  projection  factors  of  level  k  or  less  that  are  zero  in  the  cell. 
In  other  words,  if  /i, . . . ,  /^  are  the  projection  factors  of  level  k  or  less  that  are 
zero  in  c,  then  c  is  an  open  subset  of  V{<  fi, . . . ,  fs  >)■ 

For  a  (A: -I-  l)-level  projection  factor  P,  the  smallest  order  taken  on  by  P  in  c  x  M 
is  the  smallest  index  t  such  that  some  element  of  Dp,  the  set  of  all  t-th  order 
partials  of  P,  is  not  in  <  /i, . . . ,  /^  >.  The  question  of  whether  or  not  P  is  in 
fact  order  invariant  over  c  becomes  a  decision  problem: 


(3xi)  •  •  •  (3xfe)(3xfc+i) 


Fr  A 


A  [9  =  0] 

,q^Dp 


where  Fc  is  a  defining  formula  for  the  cell  c.  This  could  be  solved  by  CAD, 
or  by  other  methods.  If  the  formula  is  true,  then  P  is  not  order-invariant  over 
c,  otherwise  P  is  order  invariant  over  c.  In  fact,  when  the  above  formula  is 
true,  this  same  approach  could  be  used  to  provide  a  decomposition  of  c  x  P  into 
regions  in  which  P  is  order-invariant. 


Thus,  we  can  in  theory  deal  with  the  problem  posed  by  a  vanishing  projection 
factor.  However,  in  practice  it’s  not  clear  how  attractive  such  a  solution  would 
be.  At  the  very  least  one  may  say  that  QEPCAD  would  require  a  major  overhaul 
to  make  such  things  work. 


3.2  In  practice 

In  the  previous  section  we  saw  how  to  deal  with  the  problem  posed  by  a  vanishing 
projection  factor  ...  in  theory,  at  least.  The  question  we  consider  here  is  whether 
there  is  a  quick  and  easy  test  that  will,  in  many  practical  cases,  detect  that  the 
polynomial  P  is  in  fact  order-invariant  over  positive-dimensional  cell  c,  even 
though  it  vanishes  identically  over  c.  Another  example  will  illustrate  this  point. 
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Example  4  In  the  “edge-square  product”  problem  (see  Section  6.3),  we  have 
the  projection  factor 


P{x,  y,  Xi,X2)  =  X1X2  -  X  -1 

and  the  cell  c  =  {(— 1,?/,  0)  G  |  y  G  (/3, 7)}.  Clearly,  P  vanishes  identically 
over  c,  and  c  has  dimension  1.  However,  notice  that  y  does  not  appear  in  P. 
In  fact,  because  P  depends  only  on  x  and  xi,  and  because  the  values  for  x  and 
Xi  are  constant  within  c,  this  is  really  just  an  instance  of  our  zero-dimensional 
problem  —  Dp  =  {  — 1,  a;2,  a;i},  S  =  {— l,a;2,0},  and  the  GCD  of  the  elements 
of  S  is  1.  Thus,  P  has  invariant  order  1  in  c  x  M. 

In  practice,  cases  in  which  a  projection  factor  vanishes  identically  over  a  region 
of  dimension  greater  than  one  are  often  of  this  type.  The  hard  part  is  figuring 
out  whether  c  is  really  constant  in  all  the  coordinates  corresponding  to  variables 
that  occur  in  P. 

Our  approach  is  to  forgo  a  true  decision  procedure  for  determining  whether  c 
is  really  constant  in  all  the  coordinates  corresponding  to  variables  that  occur  in 
P  and,  instead,  settle  for  a  fast  procedure  that  is  sometimes  able  to  to  prove 
that  c  is  constant  in  all  the  coordinates  corresponding  to  variables  that  occur  in 
P.  Hopefully,  the  procedure  will  be  able  prove  that  coordinates  are  constant  in 
many  situations  in  which  this  problem  actually  occurs.  It  is  certainly  fast,  and 
is  able  to  prove  that  coordinates  are  constant  for  the  few  uncontrived  examples 
we  have  found  to  run  it  on. 

The  procedure  is  called  “CONSTCOORDTEST”.  It  takes  a  level  k,  a  fc-level 
cell  c,  and  a  polynomial  P  in  the  variables  Xi, . . .  ,Xk-  CONSTCOORDTEST 
returns  SUCCESS  if  it  is  able  to  prove  that  c  is  constant  in  all  the  coordinates 
corresponding  to  variables  that  occur  in  P,  and  FAILURE  otherwise  —  which 
means  that  it’s  unable  to  prove  this  statement  true,  not  necessarily  that  the 
statement  is  false.  Appendix  B  presents  (without  proof)  a  proposal  for  a  true 
decision  procedure  for  this  problem,  which  points  to  ways  that  CONSTCO¬ 
ORDTEST  could  be  strengthened  to  return  SUCCESS  in  more  cases. 

CONSTCOORDTEST(fc,  c,  P) 

1.  set  a  =  {}  (note:  a  grows  to  be  sort  of  a  generalized  sample  point,  where 
some  coordinates  are  given  values  while  others  remain  free,  like  (x  = 
-^,y  =  y,z  =  0)) 

2.  for  i  from  1  to  A:  do 

(a)  set  Ci  to  the  t-level  ancestor  of  c 

(b)  if  Ci  is  a  sector  then 

i.  if  degxi{P)  =  0,  add  Xi  =  Xi  to  a,  otherwise  return  FAILURE 

(c)  if  Ci  is  a  section  then 
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i.  set  L  to  the  set  of  t-level  projection  factors  of  which  Ci  is  a  sec¬ 
tion  (by  definition  of  “section”  none  of  these  polynomials  vanish 
identically  on  Ci-i) 

ii.  set  L'  to  the  set  of  all  elements  of  L  evaluated  at  Xi  =  ai,  where 
ai  is  the  ith  coordinate  of  the  sample  point  for  Ci 

iii.  set  L*  to  the  set  of  all  elements  of  L'  “evaluated  at  a”  (this  is 
actually  a  partial  evaluation,  since  some  of  the  entries  of  a  are 
of  the  form  xj  =  xj) 

iv.  if  any  element  of  L*  is  the  zero  polynomial  then  add  Xi  =  ai  to 
a,  otherwise 

if  degxiiP)  =  0  then  add  Xi  =  Xi  to  a,  otherwise  return  FAIL¬ 
URE 

3.  return  SUCCESS 


To  prove  that  this  algorithm  is  correct,  we  will  show  that  after  i  iterations,  if 
the  algorithm  has  not  returned  FAILURE,  all  the  coordinates  that  are  assigned 
values  in  a  are  “correct”  for  Cj,  in  the  sense  that  Ci  is  constant  in  those  co¬ 
ordinates  with  those  values  (although,  of  course,  ct  may  be  constant  in  other 
coordinates  as  well). 

The  i  =  0  case  is  trivial  since  a  =  0. 

Suppose  f  >  0,  consider  the  beginning  of  the  ith  iteration  of  the  loop  Step  2. 
By  induction,  a  is  “correct”  for  Ci_i  in  the  above  sense. 


Ci  is  a  sector  :  Step  2.b  determines  the  course  of  this  iteration.  If  Xi  appears 
in  P,  FAILURE  is  returned  (as  it  clearly  should  be).  If  Xi  does  not  appear 
in  P,  then  we  leave  a  unaltered  and,  because  Ci  inherits  its  coordinates 
for  {xi, . . .  ,Xi_i}  from  Ci_i,  a  is  “correct”  for  Ci  after  the  tth  iteration. 

c  is  a  section  :  Step  2.c(i)  sets  L  to  the  set  of  i-level  projection  factors  of 
which  Ci  is  a  section.  If  /  is  a  continuous  real- valued  function  defined  over 
Ci-i  whose  graph  over  Ci-i  intersects  Ci  in  at  least  one  point,  then  by  the 
usual  definitions  for  CAD: 

Ci  is  the  graph  of  /  over  a-i  {3p  G  L)(Va  G  Ci-i)[p{a,  /(a))  =  0] 

(1) 

Now,  if  f{xi, . . .  ,Xi-i)  =  (3,  and  (3  is  the  zth  coordinate  of  at  least  one 
point  in  Cj,  (1)  tells  us 

Ci  is  the  graph  oi  f  =  (3  over  Ci_i  {3p  G  L){ya  G  Ci-i)[p{a,  (3)  =  0] 

(2) 

Of  course,  “c^  is  the  graph  oi  f  =  (3  over  Ci_i”  means  that  Ci  has  constant 
ith  coordinate  with  value  l3.  Therefore,  we  have  the  following  result: 


Theorem  1  Let  (3  is  the  ith  coordinate  of  at  least  one  point  in  Ci.  Cell  Ci 
has  constant  ith  coordinate  if  and  only  if 


(3p  G  T)(Va  G  Ci_i)[p(a,/3)  =  0]  (3) 

Now,  our  “/3”  will  be  the  zth  coordinate  of  the  sample  point  for  c^. 
Thus,  Theorem  1  becomes:  Ci  has  constant  tth  coordinate  if  and  only  if 

(3p  G  L)(Va  G  Cj_i)[p(a,  Oi)  =  0] 

Step  2.c(ii)  sets  L'  to  be  the  set  of  all  p  G  T  evaluated  at  Xi  =  ai,  so  we 
may  restate  the  above  as:  Ci  has  constant  zth  coordinate  if  and  only  if 

(3(7  G  L'){'ia  G  Cj_i)[g(a)  =  0] 

Recall  that  a  is  a  “partial  sample  point”  in  which  variables  are  only  given 
values  if  the  associated  coordinates  in  Ci_i  have  been  proven  by  previous 
loop  iterations  to  be  constant.  If  we  let  “|o”  denote  “evaluation  at  a”  in 
this  sense,  then  q\a{a)  =  <7(5')  for  any  a  G  Ci_i.  Step  2.c(iii)  computes 
=  {d\a  such  that  <7  G  L'}.  So  the  above  may  be  restated  as:  ci  has 
constant  zth  coordinate  if  and  only  if 

(3g  G  L*)(Va  G  Ci_i)[(7(a)  =  0]  (4) 

Now  consider  the  test  in  Step  2.c(iv).  If  there  is  an  element  q  in  L*  that 
is  the  zero  polynomial  then  clearly  (Va  G  Ci-i)[q(a)  =  0],  which  means 
that  (4)  is  true,  which  in  turn  means  that  Cj  has  constant  zth  coordinate 
(with  value  ai).  Therefore,  if  Step  2.c(iv)  determines  that  the  coordinate 
is  constant,  it  is  correct.  Otherwise,  it  either  returns  FAILURE  or  Xi  does 
not  appear  in  P  and  it  proceeds  without  assuming  that  the  a;i-coordinate 
is  constant  in  Ci.  Thus,  after  i  iterations,  a  is  “correct”  for  Ci. 

When  CONSTCOORDTEST  returns  success,  we  find  a  minimal  delineat¬ 
ing  polynomial  for  P  just  as  we  did  in  Section  2.1  for  the  zero-dimensional 
case.  If  this  minimal  delineating  polynomial  is  constant,  we  know  that  P 
is  in  fact  order-invariant  in  c  x  M.  Thus,  once  again,  there  is  no  need  for 
QEPCAD  to  report  failure. 


4  When  is  order- invariance  really  needed 


McCallum  notes  that  when  our  goal  is  to  construct  a  sign-invariant  CAD,  as  it  is 
for  quantifier  elimination,  order-invariance  of  projection  factors  is  not  required 
at  the  highest  level.  At  lower  levels  we  need  order-invariance  to  ensure  that 
subsequent  lifting  steps  are  valid,  but  at  the  highest  level  there  are  no  subsequent 
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lifting  steps,  and  thus  sign-invariance  is  sufficient.  A  projection  factor  that  is 
identically  zero  in  some  region  c  x  K  is  certainly  sign-invariant  in  the  region,  and 
therefore  we  don’t  need  to  worry  about  delineating  polynomials  in  this  case. 

In  fact  whenever  the  projection  factor  in  question  does  not  have  a  derivation 
as  a  resultant  or  discriminant  of  other  projection  factors  only  sign-invariance 
is  required.  In  [3]  it  is  pointed  out  that  only  sign-invariance  is  required  for 
coefficients,  and  clearly  only  sign-invariance  is  required  of  input  polynomials, 
regardless  of  their  level.  In  practice  many  cases  in  which  projection  factors 
vanish  identically  involve  polynomials  that  have  no  derivations  as  resultants  or 
discriminants.  This  test  is  an  important  tool  for  recognizing  situations  in  which 
the  McCallum  projection/lifting  method  seems  to  fail,  but  in  fact  does  not. 


5  Conclusion 


The  McCallum  Projection  produces  much  smaller  projection  factor  sets  than 
previous  projections,  which  makes  it  feasible  to  attack  larger  problems  than 
would  otherwise  have  been  possible.  Unfortunately,  CAD  construction  is  compli¬ 
cated  by  the  requirement  of  sign-invariance  rather  than  simply  order-invariance 
of  projection  factors  in  cells.  This  requirement  sometimes  causes  the  McCal¬ 
lum  projection  to  “fail”,  and  sometimes  requires  the  addition  of  a  “delineating 
polynomial”  during  lifting. 

In  theory,  adding  delineating  polynomials  does  not  add  significantly  to  compu¬ 
tational  costs.  In  practice,  however,  it  complicates  implementation.  In  the  case 
of  QEPCAD,  grafting  the  facilities  for  adding  delineating  polynomials  would 
be  a  significant  undertaking.  When  the  McCallum  projection  fails,  we  have  a 
bigger  problem.  Either  another  projection  must  be  used,  which  will  likely  result 
in  a  much  larger  projection  factor  set,  or  human  intervention  will  be  needed  to 
break  the  problem  up  into  subproblems  for  which  the  McCallum  projection  can 
be  used. 

This  paper  refines  the  criteria  for  determining  when  delineating  polynomials 
are  needed,  what  delineating  polynomial  is  needed,  and  when  the  McCallum 
projection  fails.  The  result  is  that  we  can  often  safely  use  the  original  lifting 
method  with  the  McCallum  projection.  Our  implementation  in  QEPCAD  of 
these  criteria  means  that,  for  the  first  time,  the  McCallum  projection  can  be 
used  and,  in  the  absence  of  error  messages,  is  proved  to  be  valid.  Some  future 
implementation  of  CAD  will  doubtless  allow  for  the  addition  of  delineating 
polynomials  (or,  in  the  terminology  of  [3],  allow  for  adding  points  to  CAD’s). 
Such  an  implementation  would  still  benefit  from  the  criteria  developed  in  this 
paper. 
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6  Appendix  A:  Problems 


This  appendix  describes  some  of  the  basic  quantifier  elimination  problems  con¬ 
sidered  in  this  paper. 


6.1  The  x-axis  ellipse  problem 


The  x-axis  ellipse  problem,  a  special  case  of  the  general  ellipse  problem  posed 
by  Kahan  [10],  is  a  traditional  benchmark  problem  for  quantifier  elimination 
algorithms.  (See  for  example  [9],  [8])  The  problem  asks  when  the  ellipse  (x  — 
c)^  jc?  +  I =  1  lies  in  the  unit  circle.  Of  course  we  require  a  and  b  to  be 
non-zero,  and  in  fact  we  are  only  interested  in  the  case  where  they  are  positive. 
The  formula 


(Va:)(V?/) 


a>0  A  &>0  A  [b'^{x  —  cY  + 
a^b"^  =  0  — >  -I-  —  1  <  O] 


expresses  this  as  a  quantifier  elimination  problem. 


6.2  Anai’s  problem 


In  [1],  Hirokazu  Anai  applies  the  Virtual  Term  Substitution  method  of  quanti¬ 
fier  elimination  [12]  to  problems  in  control  theory.  The  particular  problem  we 
examine  is  the  last  example  posed  in  the  paper,  which  finally  boils  down  to  the 
following  quantifier  elimination  problem: 


{Exi){Ex2) 


s>0Aa;2>0Aa;i-|-l>0A  SX2  —  >  0 

Aa;2a:i  -I-  0:2  >  0  A  s(xi  -I-  1)  >  0 

A(a;2a:i  -I-  a;2)s  —  xf  —  xf  >  0  A  z  —  xi  —  X2  >  0 

A-10<sAs<0 


6.3  The  edge-square  product  problem 

Consider  the  complex  segment  L  =  {x  +  i  \  x  G  [0,  2]},  and  the  complex  square 
S  =  {x  +  iy  \  X  G  [2,4],?/  e  [—1,1]}.  Quantifier  elimination  can  be  used  to 
determine  the  complex  product  of  S  and  L.  Since  there  are  the  easily  derived 
necessary  conditions  that  the  product  lies  in  the  box  [—1,  9]  x  i[— 6, 6],  the  prod- 
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uct  can  be  expressed  as  all  pairs  (x,  y)  satisfying: 


A 


This  problem  appears  in  [4]. 


7  Appendix  B:  A  true  decision  procedure  for 
CONSTCOORDTEST 

Here  is  a  simple  algorithm,  “CONSTCOORDDECIDE”,  which  takes  a  fc-level 
cell  c  and  a,  {k  +  l)-level  projection  factor  P  and  determines  whether  or  not 
the  coordinates  of  the  points  in  c  are  constant  for  all  coordinates  corresponding 
to  variables  appearing  in  P.  It  returns  SUCCESS  when  this  is  the  case,  and 
FAILURE  otherwise.  The  algorithm  is  impractically  slow,  but  a  fast,  though 
weaker,  variant  of  it  is  the  previously  discussed  CONSTCOORDTEST. 

CONSTCOORDDECIDE(A:,  c,  P) 

1.  set  a  =  {}  (note:  a  grows  to  be  sort  of  a  generalized  sample  point,  where 
some  coordinates  are  given  values  while  others  remain  free,  like  (x  = 
-i,y  =  y,z  =  0)) 

2.  set  F  =  {} 

3.  for  i  from  1  to  A:  do 

(a)  set  Ci  to  the  Ulevel  ancestor  of  c 

(b)  if  Ci  has  dimension  0,  add  xt  =  on  to  a,  where  ai  is  the  tth  coordinate 
of  the  sample  point  for  Ci 

(c)  if  Ci  is  a  sector  then 

i.  if  degxi{P)  =  0,  add  Xi  =  Xi  to  a,  otherwise  return  FAILURE 

(d)  if  Ci  is  a  section  cell  of  dimension  greater  than  zero  then 

i.  set  L  to  the  set  of  Ulevel  projection  factors  of  which  Ci  is  a  section 

ii.  set  L'  to  the  set  of  all  elements  of  L  “evaluated  at  a”  (this  is 
actually  a  partial  evaluation,  since  some  of  the  entries  of  a  are 
of  the  form  Xj  =  Xj),  removing  any  that  are  identically  zero 

iii.  set  L*  to  the  set  of  all  elements  of  L'  evaluated  at  Xi  =  ai,  where 
ai  is  the  ith  coordinate  of  the  sample  point  for  Ci 


(3xi)(3x2)(3y2) 


X  =  xiX2  -y2  A  y  =  xiy2  +  X2 
0<xi<2  A  2<a;2<4A 
-1  <  2/2  <  1 

A  — l<a:<9  A  — 6<2/<6 
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iv.  if  any  element  of  L*  is  in  the  ideal  defined  by  the  irreducible  com¬ 
ponent  of  V (F)  that  contains  Ci  (note  that  either  all  elements  of 
L*  are  in  the  radical  of  the  ideal  generated  by  the  elements  of 
F,  or  none  are)  then  add  Xi  =  ai  to  a,  otherwise 

if  degj;.{P)  =  0  then  set  F  =  F  LI  L'  and  add  Xi  =  Xi  to  a, 
otherwise  return  FAILURE 

4.  return  SUCCESS 


As  a  simple  example,  let  P{xi,X2,  X3,X4,  z)  =  {xi  —  X4  +  l)z‘^  —  xi  and  let  c  be 
the  1-dimensional  cell  in  constructed  via  the  following  lifting  steps: 


cell  type 

sections  of  ... 

sample  point  coord. 

Level  1 

section 

{a;i} 

oi  =  0 

Level  2 

sector 

— 

02  =  3/2 

Level  3 

section 

{2:3  -1-  X2} 

03  =  -3/2 

Level  4 

section 

{2:3X4  +  X2+  Xi} 

04  =  1 

At  the  start  of  the  4th  iteration  of  CONSTCOORDDECIDE,  F  =  {0:3  -I-  X2} 
and  a  =  {xi  =  0,X2  =  X2,X3  =  X3}.  First  L  is  set  to  {X3X4  +  X2  +  Xi\,  then 
L'  is  set  to  {X3X4  +  X2},  and  finally  L*  is  set  to  {0:3  +  X2}-  Clearly,  the  one 
and  only  element  of  L*  is  in  the  ideal  defined  by  the  one  and  only  component 
of  U(F),  and  therefore  coordinate  4  of  every  point  in  c  has  the  value  1. 

In  practice  CONSTCOORDTEST  seems  to  be  reasonably  successful  at  recog¬ 
nizing  when  a  cell’s  coordinates  are  constant  at  each  coordinate  corresponding 
to  a  variable  that  appears  in  P.  This  can  perhaps  be  explained  as  follows:  The 
elements  of  F  are  polynomials  in  the  variables  that  do  not  appear  in  P.  The 
discriminant  of  P  and  other  projection  factors  descended  from  P  do  not  con¬ 
tain  these  variables,  so  if  they  appear  in  L,  then  the  ideal  membership  from 
Step  3.d(iv)  reduces  to  simple  zero  testing,  which  is  precisely  what  CONSTCO¬ 
ORDTEST  does. 
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